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Three-dimensional quantum phase diagram of the exact ground states of a mixture of 
two species of spin-1 Bose gases with interspecies spin exchange 
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We find nearly all the exact ground states of a mixture of two species of spin-1 atoms with 
both interspecies and intraspecies spin exchanges in absence of a magnetic field. The quantum 
phase diagram in the three-dimensional parameter space and its two-dimensional cross sections 
are described. The boundaries where the ground states are either continuous or discontinuous are 
determined, with the latter identified as where quantum phase transitions take place. The two 
species are always disentangled if the interspecies spin coupling is ferromagnetic or zero. Quantum 
phase transitions occur when the interspecies spin coupling varies between antiferromagtic and zero 
or ferromagnetic while the two intraspecies spin couplings both remain ferromagnetic. On the other 
hand, by tuning the interspecies spin coupling from zero to antiferromagnetic and then back to zero, 
one can circumvent the quantum phase transition due to sign change of the intraspecies spin coupling 
of a single species, which is spin-decoupled with the other species with ferromagnetic intraspecies 
spin coupling. Overall speaking, interplay among interspecies and two intraspecies spin exchanges 
significantly enriches quantum phases of spinor atomic gases. 

PACS numbers: 03.75.Mn, 03.75.Gg 



I. INTRODUCTION 

Spinor Bose gases have been a remarkable subject since 
a decade ago . A mixture of two scalar Bose gases 
has also been studied extensively @, ^ . The problem of 
mixing two distinguishable species of spinor Bose gases 
with interspecies spin exchange was first considered in 
a mixture of pseudospin-i gases, with particular mo- 
tivation and interest in entanglement between the two 
species, realizing the so-called entangled Bose-Einstein 
condensation [l3-[l3|, which is a feature beyond those of 
scalar mixtures as well as a single species of spinor gases. 
A mixture of two species of spin-1 atoms has also re- 
ceived attention more recently [l3 - [l7l |. Compared with 
a pseudospin-i mixture, the spin-1 mixture may be more 
experimentally accessible and exhibits interesting prop- 
erties such as the coexistence of intraspecies and inter- 
species singlet pairs even in a global many-body singlet 
state. Previously, the exact ground states of a spin-1 mix- 
ture have been found in four parameter regimes, with or 
without a magnetic field [16j]. This left open to work 
out the ground states in other parameter regimes, espe- 
cially, an overall picture of the occupation of different 
ground states in the whole parameter space. However, 
the present problem of minimizing the energy is very 
complicated, with three constrained variables, namely, 
the total spins of the individual species and of the total 
mixture, as well as four independent parameters, namely, 
the three spin coupling strengths and the magnetic field. 

In this paper, we find nearly all the possible ground 
states in absence of a magnetic field, by working out the 
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complicated minimization problem (the even more dif- 
ficult case with a magnetic field is discussed in a sepa- 
rate paper). Thus we obtain a three dimensional phase 
diagram, with the three spin coupling strengths as the 
parameters. By comparing the ground states in neigh- 
boring parameter regimes, we distinguish phase bound- 
aries, where discontinuities of the ground states or quan- 
tum phase transitions occur, from crossover boundaries, 
where the ground states are continuously connected. 
As the spin coupling strengths, proportional to spin- 
exchange scattering lengths, can be tuned through Fesh- 
bach resonances, a novel venue is opened in this system 
for studying quantum phase transitions and crossovers 
between different many-body ground states. 

Consider two species a and b of spin-1 atoms, whose 
numbers Na and Ni, are conserved respectively. The 
many-body Hamiltonian has been given previously [l6| . 
In absence of a magnetic field, it can be written as 

^ = 1 S^ + ^S^ + c-'-S^.S,, (1) 

where — ajj^Ff^^a^, is the total spin operator for 
a species, is the annihilation operator of species a 
(a = a, 6), F^^ is the {^J.l') element of spin-1 matrix 
— —1,0,1), c" is the intraspecies spin coupling 
strength of species a, proportional to the intraspecies 
spin-exchange scattering length, while c"''' is the inter- 
species spin coupling strength, proportional to the inter- 
species spin-exchange scattering length [iGj . 

Sa , Sb and the total spin S are all good quantum num- 
bers. Consequently, there are 25* -I- 1 degenerate ground 
states with = -S"™, • • • , 5™, i.e. 

I rrm nm Qm q \ 

where 5™, 5™ and 5™ are, respectively, the values of Sa, 
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Sb and S that minimize the energy 

a _ ab b _ ab 
E = SaiSa + 1) + SbiSb + 1) 

2 .b 2 (2) 

+ ^5(5 + 1), 

under the constraints 

< 5a < iVa, (3) 
|5a-5b| < 5< 5a+56. (4) 

In general, the problem of minimizing E is very com- 
plicated, because there are three variables Sa, St and 5, 
with the constraints ([3|) and ([4]), as well as three param- 
eters. Through tour de force calculations, we have done 
the minimization in nearly all parameter regimes, except 
a small regime where there are too many possibilities. 
The details of solving for 5'™, 5™ and S*™ are reported 
in the Appendices. 

In discussing the phase diagrams, we shall focus on the 
maximally polarized state 

5 m rtm cm nrnx 
, a ' '-'6 ' 1 '-' /i 

which can be determined to be the unique ground state 
under an infinitesimal magnetic field. 

In Sec. ini we list the ground states in different pa- 
rameter regimes in a table and draw the phase diagrams 
on c" — planes for different signs of c"'', based on 
calculations described in the Appendices. In Sec IIIIl we 
describe the structure of the phase diagram in the three- 
dimensional parameter space, and its implication on ef- 
fects of c"'' on quantum phase transitions. In Sec. IIV[ 
we discuss composition of the ground states in terms of 
bosonic degrees of freedom. A summary is made in SecfVl 

II. GROUND STATES IN VARIOUS 
PARAMETER REGIMES 

The ground states in different parameter regimes with 
c"* ^ are reported in Table HI Different parameter 
regimes are numbered according to the ordering numbers 
of the sections and subsections in the Appendices where 
the minimization of E is done in these different regimes 
respectively. 

Although the calculations are very complicated, the re- 
sults are quite elegant, as displayed in the phase diagrams 
shown below. 

A. c"'' < 

For ferromagnetic interspecies spin coupling c"'^ < 0, 
the maximally polarized ground state is in the form of 

I cm cm cm i om cm i c'^^\ I c^\ I Q^^ c^" \ 

Pa T '^b T '^a + '^b T^a Pa i '-'a /aPb i ^b /b, 

(5) 



which is always disentangled. However, the disentangled 
ground state may still be different in different parameter 
regimes. We can determine whether the ground states 
in neighboring regimes belong to a same quantum phase 
from whether they are continuously connected, that is, 
whether they both approach the ground state on the 
boundary between their regimes. 

In FIG.[l] we draw the two-dimensional phase diagram 
of the maximally polarized ground state for c"*" < 0. 

The ground state is |0, 0)o|0, 0)f, in the first quadrant 
c'^ > and > excluding the part < < 

and < c'' < '^wt^- 

In the regime c" < =^^^ and c'' < '^N.ti" ' 
eluding the whole third quadrant, the ground state is 

\Na,Na)a\Nb,Nb)b. 

In the regime Ca < and < Cf, < -2A^Qc"^ the 

ground state is \Na, Na)a\n'i,ni)b, which varies with c'', 

as ni = Int[ ^'^ — ^] depends on c'^'' and c''. It is con- 
tinuously connected upwards with \Na, Na)i\0,0)b on the 
boundary c'' — —2Nac"-'', and is also continuously con- 
nected downwards with \Na, Na)a\Nb, Nb)b on the bound- 
ary c*" = - ^7 ■ Hencetheregimeof |7Va,7Va)a|ni,ni)f, 
acts as a crossover regime. 

Likewise, in the regime < and — 2n'''+i — < 
~2NbC°-^^ the ground state is 1^2, n2)a| A'^f,, A^f,)^, which 
varies with c°, as n2 = Int[ ^'^ ^J^* — i] depends on 
and c°. Throughout the paper, Int(a;) refers to 
the integer closest to x and in the legitimate range 
of the concerned quantity. This regime is continu- 
ously connected leftwards with \NaTNa)a\Nb,Nb)b on 
the boundary = ~ 2jv''+i ' ^^"^ ^'^ continuously con- 
nected rightwards with |0, 0)a|A^b, Af;,)f, on the boundary 
= -2NbC°-^. Hence the regime of |n2, n2)Q|A^6, A^6)& 
acts as a crossover regime. 

Altogether, there are only two phases for c"'' < 0. 
One includes these five ground states in the second, third 
and fourth quadrants, including the borders = and 

= 0, plus the regime < c° < and < c'' < 

— I^q-j-. The other is just the ground state |0,0)Q|0,0)b 
in the complementary regime. A discontinuity or quan- 
tum phase transition occurs in crossing anywhere on the 
boundary between the two phases, which is a continuous 
line consisting of the intervals > while c" = 0, 

< < c"^ while & = m^^, < < c"'' while 

c° = i^TT' ^'^ ^ i^TT ^^^""^ = 0. The ground 
state of the total system being a product of the states of 
the the two species, the state of one or both species dis- 
continue on this phase boundary. These quantum phase 
transitions are first order, as the total Sz is discontinuous 
in crossing anywhere on the the phase boundary. 

The ground states on the phase boundary itself need 
to be addressed in detail. In the present case of c"'' < 0, 
we have carefully taken into consideration the boundaries 
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No. 


Parameter regimes 


Ground st&tcs 


l,A2a 


c"'' < 


c" < 

b ^ 2NaC°-'' 
- 2JVh + l 


\Na,Nb,Na + Nb,S,), 

disentangled when Sz = ±{Na + Nb), entangled otherwise 


A2b 


c" < 0, 

<Cb< -2N,c'^» 


\N,,nuN,+nuSz), m = Int[J^^^ - i], 

disentangled when Sz — ±{Na +ni), entangled otherwise 


A2c 


< 0. 

<!• > -2A-„r"'' 


\N„.O.K,,S,) = \N„,S,)„\0.0}b, 

l.Lii)',. ill, (lli^illj. Vvliv. ii 1. ' —I : ' (J • "---ii I f lii J-,iv t. I V ' I \\ iril,^ 


l,A3a 


c" < 0, 

- 2JV„ + 1 


lAa,Af„Aa. + Ai^S^), 

disentangled when Sz = i{Na + Nb), entangled otherwise 


A3b 


c'' < 0, 


\n2,Nb,n2 + Nb,Sz), n2 = Int[l^;^ - i], 
disentangled when Sz = ±(n2 + Nb), entangled otherwise 


A3c 


c" < 0, 

c" > -2Nbc'''' 


\0,Nb,Nb,Sz), 

disentangled when Sz = iNb, entangled otherwise 


A4 


c" > 0, < c" < -2iV„c"^ 


|0, 0,0,0) = |0,0)a»|0,0)6, 
disentangled 


A5 


c" > 0, < c" < -2iV^,c"^ 

c c > [c ) 


A6 


C > U, C > —ZNaC 


A7 


c" > 0, c° > -2Ai,c'"' 


A8 


c" > 0, -TKfxT < c*" < -2A^„c'''' 


A9 


> ^' 2]V„ + 1 ^ C < 2iV6C , 


AlO 


o<c < 2jv;+i' 
< c < 2iv;+i 


\Na,Nb,Na+Nb,Sz}, 

disentangled when Sz = ±(n2 + Nb), entangled otherwise 


Bl 


c"'' > 0, 

Na > Nb 


a b ab 
C = c = C , 


\sr,sr,o,Q),sr = o,---,Nb 

always entangled 


B2a 


C < C , 
„a ^ 2JVi,c»'' 
— 2JV„+1 


\Na,Nb,Na-Nb,Sz), 
dilWciyo cIltcLll^itJU. 


B2b 


2]Va + l — ^ — 2Ar^+l ' 


\n2,Nb,n2 - Nb,Sz), 
always entangled 


B2c 


c ^ c , 

aiVfcc"'' ^ a ^ ab 
2JVft + l — — ' 


\Nb,Nb,0,Q) 
entangled 


B345 


C„ > 0. > 0, <:"'' > 


|n. 0,0.0) = |o.o)„ « |o.o)6 

(lis(.'iitaiigl(.'( I 


B6a 


< 0, 

„ab ^ a ^ 2(]Vt + l)c°'' 
^ ^ — 2JVh + l 


|Af„A5,0,0) 
entangled 


B6b 


^b ^ r, 

^ £ 2(A, + !)<:'"■ 


Ins, Nb, Nb - ns, S,), ns = Int[^^^^^^ - i] 
eiitaugk'd uuloss r" = 2(A'j; + i.e. iiy, = 


B6c 


c° < 0, 

c" > 2(Ari, + l)c"'' 


\0,Nb,Nb,Sz) = \0,0)a\Nb,Sz)b 
disent angled 


B7a 


c" < 0, 

,.ab ^ b ^ 2(JV„+l)c"'' 
^ — 2iVh+l 


\Na,Nb,Na-Nb,Sz} 

always entangled 


B7b 


c" < 0, 

'^^M^+f < < 2(A. + l)c"'' 


|A^„,n4,A^a -n4,S.), n4 = Int['^°+,^'''"" - i] 
entangled unless c*" = 2(A'a + l)c''\ i.e. n4 = 


B7c 


c" < 0, 

> 2{Na + 1)0"" 


\Na,0,Na,Sz} = \Na,Sz}a\0,0}b 

disentangled 


B8 


c" > 0, c" > 0, 
c" + c*" < 2c"'' 


\Sr, Sr, \ST - SnSz), ST ^ 0, Sr ^ 
always entangled 



TABLE I: Ground states of a mixture of two spin-1 atomic gases with interspecies spin coupling c" 7^ in different parameter 
regimes in the absence of a magnetic field. For c"*" > 0, A^a > Nb is assumed without loss of generality. Int(x) refers to the 
integer closest to x and in the legitimate range of the concerned spin quantum number, that is, < ni < Nb, Q < n2 < Na, 
0<n3<Nb,0<n4< Na. 



4 



C : 

\Na,Na)a\0.0)b 


|0,0)a|0,0),, 


\Na,Na)a\ni,ni)i 


\AT AM 1 AL AT, \ 


5 \n2,n2)a\Nb,Nb)b 


c 

|0,0)a|A^f„iVfa) 



FIG. 1: Ground states in c" — c*" parameter plane for c"*" < 0. |0,0)a|0,0)6 discontinues with all the other ground states, 
which belong to a same quantum phase, as they are continuous on the boundaries between them. For brevity, the coordinates 
of the boundaries are not labeled, as they are clearly stated in the main text and Table HI The phase boundaries belong to the 
phase other than |0, 0)alO, 0)6. 



in the calculations described in Appendix A. It turns out 
that the ground states on the phase boundary are always 
the same as the ground states on the left or downside, 
that is, in the phase other than |0, 0)a|0, 0)f,. 



B. 



For antiferromagnetic interspecies spin coupling c > 
0, the maximally polarized ground state is in the form of 

\ST,Sl\\ST~Sl%\ST~Sr\), (6) 

which is always entangled if S"^ ^ and S*™ 7^ 0. If at 
least one of them is 0, then the state is disentangled as 
|S™,Sr)a|0,0)b or |0,0)a|^r,^r)fc- In FIG. H we draw 
the two-dimensional phase diagram of the maximally po- 



larized ground state l^^", 5^", \S^ - S^\,\S^ - S^l^l) for 
c'^^ > 0. 

Without loss of generality, we assume Na > Ni, in 
the case c"'' > 0. The ground state is \Na,Nt,Na - 
Nb,Na- Nb) in the regime c'' < and < 

plus < and c"''' < < '^^I'^^l^f ■ It is continu- 
ous connected upwards with \Na, n^, Na — n^, Sz) in the 
regime < and ^^^"j^^lY" < < 2(iVa + l)c''''. 

With 71,4 = Int[ ^^ '^P'^ ^] varying with from Nb 

to 0, \Na,n4, Na — n4,Sz) is continuously connected up- 
wards with |iVa, iVa)a|0, 0)b in the regime c° < and 
c'' > 2{Na + l)c'^^ 

To the right, \Na,Nb,Na - Nb,Na - Nb) is continu- 
ously connected with \n2jNb,n2 — Nb,n2 — Nb) in the 
regime < c"^ and < c" < ^pf^- As 712 = 



i] varies from Na to Nb, depending on c", 
A^6; ^2 — Nb) is also continuously connected 



Int[ 

\n2,Nb,n2 

rightwards with | A^f,, A''b, 0, 0) in the regime 



and 



2Afi,c'' 



(excluding c° = c° = c ) 



The latter 



plus < Q and c"^ < c" < ^''^"j^^^X" 
part is further continuously connected rightwards with 
\n3, Nb, Nb — n^jNb — n^) in the regime c'' < and 

&hi)£!! <r .a ^ 9^^^, _Ln.a6 - Int[ "'"'^^' + ^^ 



< < 2iNb + l)c'^^ As ns 
i] varies from Nb to 0, \n3,Nb,Nb — n^,Nb 



n^) is also 

continuously connected rightwards with \Q,0)a\Nb,Nb)b 
in the regime c*" < and c°- > 2{Nb + l)c°''. Note 
that the apparent asymmetry between a and b is due 
to Na > Nb by definition. If Na — Nb, then the phase di- 
agram structure is indeed symmetric between a and b, as 
the regime of |7i2, Ah, 712 — Nb, n2 — Nb) disappears, while 
|A^b, Nb, 0, 0) = |A^a, A^b, Na - Nb, Na - Nb) and thus their 
regimes merge. 

Therefore, for > 0, on c° — parameter plane, 
the second, third and fourth quadrants, plus the regime 



< c"'' and c'' < c"'' (excluding c" 



) in 



the first quadrant, form a single quantum phase. Within 
this phase, an entangled state jA'a, A'f,, A'q — Nb,Na — 
Nb) on the lower left part crossovers upwards finally to 
the disentangled states IA'q, A^q)q|0, 0)^, and rightwards 
finally to the disentangled state |0, 0)a|A'6, Ah)f,. 

This regime of \Nb,Nb, 0,0) neighbors the regime 
c'^c^ > (c''^)^, c°- > and > 0, where the ground state 
is |0, 0)a|0, 0)b, on c°- — c'' — c""^ , where there are degener- 
ate ground states |S'^, S"^", 0, 0) with S"^ = 0, • • • , A^, two 
of which are the ground states in the two regimes respec- 
tively. Starting from one of these two regimes and varying 



the parameters towards the boundary c" = c 



„a6 



, the 
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\Na,Na)a\0,0)b 



\Na, n4, Na - n4, Na - Ui) 



|0,0)a|0,0)fc 



\Na,N^„Na- Nt,,Na- Ni) U^, Nb, Nb - 




\0,0)a\Nb,Nb)b 



FIG. 2: Ground states in c" — c'' parameter plane for c"'' > 0. All the ground states in the second, third and fourth quadrants 
plus the regime c" < c"'' and c*" < c"'' in the first quadrant are continuous connected, forming a same quantum phase, which 
discontinues with |0,0)a|0,0)b at = c** = c"''. The ground states in the regime above c" < c"' and c'' < c"'' and below 
c"'c'' < (c"'')^ in the first quadrant have not determined, although it is certain that those in the subregime c° + c' < 2c"'' 
in the first quadrant is always entangled. Limited by space, |n3, A^i,, A^;, — n3,Ni, — n^) is indicated by n3,Nt,Ni, — na. It is 
continuously connected leftwards with \Nb, Nt, 0, 0), which is continuously connected leftwards with |n2, Nt, n2 — Nb, 712 — A'^i,), 
which is continuously connected leftwards with \Na, Nt,, Na — Nt,,Na — Nt,}. 



ground state remains as the starting one on the bound- 
ary c"" — c'' — c""^ , and then jumps to the ground state 
in the other regime upon entering it. As S' = 5^ = in 
both |iV{,, A^b, 0, 0) and |0, 0)q|0, 0);,, the quantum phase 
transition is a continuous phase transition. 

Unfortunately, we did not work out the detailed ground 
states in the part of the first quadrant with c°c'' < (c"'')^ 
while c° > c"'' and c** > c"'', as the situation is too compli- 
cated. However, we have shown in the Appendix IB 81 that 
in the sub-regime with c°- > 0, > Q and c°- + < 2c°-^, 
the ground state is always the entangled state of the 
form \S^,Sr, |5™ - S^l \S^ - S^\), where ^ 
and 5™ ^ 0. This is in consistency with the states 
\Na, Nb, Na - Nb, Na ~ Nb) and \n2,Nb, Nb - 112, Nb - ns) 
whose regimes overlap with the present one. 

c°'S' = (c"'')^ is an equilateral hyperbola with c° and 
coordinate axes as the asymptotes, and the length of 
the real axis 

C. c"'' = 

If c"'' — 0, then % = Ha + Hb, where Ha is the a- 
species part of the Hamiltonian {a — a,b). Consequently 
the ground state is the direct product 15™, Saz)\SJ^, Sbz) 
of the respective degenerate ground states of Ha and Hb, 
with 5™ and 5™ determined by minimizing them in- 
dependently. The maximally polarized ground state is 
15"^, S'")|5^, S^) is of the same form © as in the case 
of c""^ < 0. 



As depicted in FIG. |3l the phase diagram of the maxi- 
mally polarized ground states for c"'' = simply consists 
of four quadrants, in each of which the ground state is a 
product state, that is, |0, 0)a|0, 0)^ for c° > and > 0, 
\Na,Na)a\0,0)b for c'^ < and > 0, \0,0)a\Nb, Nb)b 
for c° > and < 0, and \Na,Na)a\Nb, Nb)b for c" < 
and c** < 0. When c" = 0, the ground state of species 
a is a degenerate one \Sa, Sa), with any legitimate value 
of Sa , while the ground state of the other species is still 
determined by the sign of its intraspecies spin coupling. 
With c"-^ — 0, the origin c° = c'' = is a quadra-critical 
point, where the ground state is any state of the form 
of I Sa, Sa)\Sb, Sb) , degenerate in both Sa and Sb, or any 
of the superposition of these states. With the disconti- 
nuities of ground state and of the total Sz , a first order 
quantum phase transition occurs in crossing anyway on 
the boundaries. If the ground state is not limited to 
be maximally polarized, then it can be \Sa, Saz)\Sb, Sbz) 
with arbitrary legitimate values of Sa , Saz , Sb and Sbz or 
any of their superposition. 



III. THREE DIMENSIONAL PHASE DIAGRAM 

A. Structure 

With the two-dimensional c° — c'' phase diagrams 
given above for all possible values of c"'', we actually have 
obtained the three dimensional phase diagram with c". 
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\Na,Na)a\0,0)b 


|0,0)a|0,0)fc 


\Na,Na)a\Nb,Nb)b 


( 

\0,0)a\Nb,Ni,)b 



FIG. 3: Ground states in c" — c' parameter plane for c"'' = 0, which are direct products of the ground states of the two species. 
If c" = 0, then the ground state of species a is degenerate with any legitimate value of 5*^, (a — a, b). 



and c"'' being the three coordinates. It is very difficult 
to draw such a three-dimensional phase diagram, so we 
only make some descriptions in the following. 

In the three dimensional parameter space, a boundary 



= "fc""^ represents a plane intersecting the plane c"'' = 
on c" axis and with an angle arctan(l/7). Here a ^ fi 
represent the two species, 7 represents a coefficient. 

In the half space of c"'' < 0, the phase boundary sepa- 
rating |0, 0)a|0, 0);, from the other phase is given by the 

planes Cb = — 2jvt+i ^"^"^ ~ ~ 2jv''+i ^^^^ octant 

(c° > 0, c'' > 0, c°-^ < 0), together with the plane c"" = 

for Cb > - ajVt+i ^'^'^ t^*^ plane c'' = for Ca > -f^^- 
Within the other phase, the crossover boundaries are 

planes Cb — — 2jvt+i ^^'^ ~ — 2A^aC°^ in the sixth oc- 
tant (cfe > 0, Ca < 0, c"'' < 0), and planes Ca = - 2jvl+i 
and Ca = —2NbC°'^ in the eighth octant (cq < 0, Cb < 0, 
c'^'' < 0). 

Now consider the half space c"'' > 0. The bound- 
ary surface of |0, 0)a|0,0)b is a quarter of the hyperbolic 
paraboloid c°'c^ — (c"'')^, with all the three coordinates 
positive, i.e. in the first octant. Its intersection figure 
with a fixed value of gives the hyperbola boundary 
in the two-dimensional diagram as discussed in the last 
section. The intersection figure with a fixed value of c** 
is half of a parabola with the vertex at c" = c"'' = 0, i.e. 
on the c^ axis, and focus distance c''/4. The intersection 
figure with a fixed value of c"" is half of a parabola with 



the vertex at — c°-^ 
distance c"/4. 

For c°-^ > 0, the crossover boundaries c° = 2n'''+i ' 

M>>^ and c° = c"'' under the constraint < c"'', 

and c" = 2{Nb + l)c°'' under 



^ 2Ni + l 

as well as c° 



0. i.e. on the c° axis, and focus 



the constraint c'' < are all planes starting from, but 
excluding, c^ axis and extending to infinity. Likewise, 



the crossover boundaries c — 



J, _ 2(Ar„-H)c° 



2Nf, + l ' 

2(A^a + l)c°'' under the constraint c° < are all planes 
starting from, but excluding, the c°' axis and extending 
to infinity. 

c° = c'' = > 0, where the two phases neighbor and 
there are degenerate ground states 15^,5^,0,0), now 
represents a straight line starting from the origin and 
extending in the first octant. c° -I- c** < 2c°^, > and 

> now represent the part of the first octant that is 
surrounded by plane c"" + c^ — 2c"'', which starts from 
the origin, together the planes c° = and c^ — 0. 

In FIG. m we draw the two-dimensional c^ — c"'' cross 
section for a given < 0. The intersection figures 
with the continuous connecting boundaries are the four 
oblique lines starting from, but excluding, <f ~ c"^'' = 
(i.e. c° axis). The ground states on the two sides of each 
of these four boundaries arc continuously connected. But 
\Na,Na)a\Nb,Nb)b and \Na,Nb,Na - Nb,Na - Nb) dis- 



continue on the negative half of c axis (c < 0, 



0). 



2Nb + l 



At c = c = 0, the maximally polarized ground states 
are degenerate, being \Na, Na)a\Sb, 5f,)h, where Sb is any 
legitimate spin quantum number of b species. 

Similarly, in FIG. [HI we draw the two-dimensional 
c° — c""^ cross section for a given c*" < 0. The intersec- 
tion figures with the crossover boundaries are given as 
six oblique lines (excluding the origin c'' — c""^ — 0). The 
ground stats on the two sides of each of these six bound- 
aries are continuously connected. But \Na,Na)a\Nb,Nb)b 
and \Na,Nb^ Na — Nb, Na — Nb) discontinue on the neg- 
ative half of c" axis {c" < 0, c"'' = 0). At = c"'' = 0, 
the maximally polarized ground states are degenerate, 
being \Sa, Sa)a\Nb, Nb)a, where Sa is any legitimate spin 
quantum number of a species. 
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\Na,Na)a\ni,ni)^ 



c}" = 2{Na + 1)0'^'' 



\Na,n4,Na - n4,^a " "■4) 




\Na, Na)a\Nb, Nb)b l^a, Nb, Na " iVfe, Na " iVfc) 



FIG. 4: c'' - c"^ phase dia gram for a given c" < 0. The ground state is continuous on the four crossover boundaries represented 
as the four obhque lines (excluding c'' — c"'' = 0), but discontinues on the negative half of c*" axis. 



|n2,n2)a|^fe,-^Vb)&^ 




n _ 2NbC' 



|"2, Nb, n2 - Nb, n2 - Nb) 

,a ^ 2(jVt, + l)c°'' 

\Nb,NbM 



|n3, Nb, na - Nb, - Nb) 



Na, Na)a\Nb, Nb)b \Na, Nb, Na " Nb, Na - Nb) 



FIG. 5: c" - c"' phase dia gram for a given c** < 0. The ground state is continuous on the four crossover boundaries represented 
by the six oblique lines (excluding c" — c"'' = 0), but discontinues on the negative half of c*" axis. 



B. Convergence of the boundaries as c" 



ent ground states converge to those of c = 



In the two-dimensional — c phase diagram with each 
value of 0°"^ , the distance from c° = c'' = to each bound- 
ary between two ground states are proportional to 

In the three-dimensional parameter space, as d 



^ah I 



.ah 



0, 



However, as 



^ah 



0, the crossover regimes in c° 



all boundaries evolve to c" — c axes. All boundaries c phase diagrams diminish, and the same phase in the 

second, third and fourth quadrants tend to be separated 
to three phases. Indeed, as we know, there is no crossover 
regime in c° — phase diagram for c"'' = 0. In a c° — 



of the form c" = evolve toward c" = c = 0, i.e. 
axis (/3 ^ a), d^c^ — (c"'')^ evolves towards and 
axes. The lines c" ^ ^ c"'' and + = 2d-^ in 

the first octant both evolve towards the origin c" = c*" 

d 



= 0. 



cross-section for < (FIG. H]) and in a 



^ah 



Therefore, as d 



.ah 



0, the boundaries between differ- 



cross-section for d < Q (FIG. O, the boundary lines of 
continuous connection converge to the respective origins. 
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C. Quantum phase transitions caused by varying 

ab 

C 

The ground states on = plane is always contin- 
uous with those for c"'' < 0, because when c"'' = —e, 
where e represents a positive infinitesimal quantity, the 
ground states in each regime of c° and is the same as 
^ab _ Q Moreover, in first, second and fourth quadrants, 
the ground states on c'^^ = plane is also continuous 
with those for c"'' > 0. 

The only exception is that the ground state on the 
third quadrant of c"'' — plane is discontinuous with 
that for c°'^ > 0, although the boundary of the c" — 
regime is continuous from c"''' > to c"'' — 0. The ground 
state is \Na, N^, Na - m, Na ~ m) if Cab > 0, c'' < and 
< 0, but is \N,,N,)a\Nb,Nt)b if Cab < 0, c'' < 
and c'' < or degenerate product states if Cab and one 
or both of c° and c** are zero. Hence there is always a 
discontinuity when c"'' changes its sign from positive to 
zero or negative, while and c'' are negative or zero. 

Therefore there is a quantum phase transition in cross- 
ing anywhere on the quarter plane c"'' = 0, c" < 0, Cfc < 
(strictly speaking, between c"^ > and c"'' < 0), depicted 
as the negative half of axis in FIG. [3] and the nega- 
tive half of c'' axis in FIG. [5l This is the only regime 
where quantum phase transitions are due to variation of 
interspecies spin coupling strength. This quantum phase 
transition is first order, as there is a discontinuity of Sz, 
which is Na + Nb for c'^^ < 0, but is \Na-Nb\ for c"'' > 0. 

However, in the three-dimensional parameter space, 
these two ground states can be continuously related 
through the crossover regimes with c° < and c*" > 0, 
as indicated in FIG. IH or through the crossover regimes 
with c'' < and c° > 0, as indicated in FIG. [3 In 
other words, the quantum phase transition can be cir- 
cumvented if the sign of c"''' is changed while one of c° 
and is changed from negative to positive and then back. 

D. Avoiding single-species quantum phase 
transition by tuning c"'' 

In the three-dimensional parameter space, the part of 
the first octant with < c° < and < < c"* 
(excluding c° = c*" c"''), the part of the fifth octant 

with < < -|^^ and < c'' < -^f^, and 
all the other six octants correspond to a same quantum 
phase. The reason is the following. 

The ground state \Na,Na)a\0,0)b occupies the regime 
c'^ < and c*" > -2Nac'''' in c"'' < half space, the 
whole second quadrant of the c"^ — plane, as well as 
the regime < and c'' > 2{Na + l)c"''' in c"''' > 
half spcLce. In other words, No) a\^i^)b occupies the 
regime -c''/{2Na) < c"^ < cV[2(iVa + 1)] in the < 
and > subspace, fully including the second quadrant 
of the c"'' = plane (see FIG. H]). Similarly, the ground 
state |0, 0)a|-/V6, A^b)6 occupies the regime —c°'/{2Nb) < 
c"^ < c°-/[2{Nb + 1)] in the c" > and c'' < half space, 
fully including the fourth quadrant of the c"^ = plane 



(see FIG.©. 

Consequently some quantum phase transitions that are 
inevitable when c'^^ — can be avoided by introduc- 
ing negative c"''. With c"'' = 0, there is a first order 
quantum phase transition when c*" changes its sign while 
c-^ < 0, from \N,,Na)a\0,0)bto \Na,Na)a\Nb,Nb)b, as in- 
dicated in FIG. [3] (see also c'' axis in FIG. |4]). These two 
ground states can be continuously related through the 
crossover regime with c° < and c"'' < in the three di- 
mensional parameter space (see FIG. S]). Similarly, with 
^ab _ there is a quantum phase transition when c° 
changes its sign while < 0, from \0,0)b\Nb, Nb)a to 
\Na,Na)a\Nb,Nb)b, as indicated in FIG. [3] (see also 
axis in FIG. [5]). These two ground states can be contin- 
uously related through the crossover regime with < 
and c"'' < in the three dimensional parameter space 
(see FIG. ED. 

This implies that the quantum phase transition of a 
first species of spin-1 atoms, in varying the intraspecies 
spin-exchange scattering length between negative and 
positive, can be circumvented by mixing it with a second 
species of spin-1 atoms with negative intraspecies spin- 
exchange scattering length, and then tuning the inter- 
species spin exchange scattering length from zero to neg- 
ative and then back, in the same time as the intraspecies 
spin-exchange scattering length of the first species varies. 

IV. COMPOSITE STRUCTURES IN BOSONIC 
DEGREES OF FREEDOM 

Any state \Sa, Sb, S, Sz) of the mixture can be written 
in terms of bosonic operators, as 

\Sa, Sb, S, Sz) — ''^g{Sbz)\Sa, Sz — Sbz)a\Sb, Sbz)b, (7) 
St. 

where 

I Set : Saz} a — [r(5a)--T(^a.-fl)]-l(^a_)^°-^°^| 1 Set) a •) 

(8) 

with 

Sa_ = \/2(a|Q!o -faja-i), (9) 
r(ma) = \/ {Sa + ni){Sa -ma + 1), (10) 

and 

subject to the condition that Na — Sa is even, 

et ^ a\al, - aj' (12) 

is the creation operator of an intraspecies singlet pair of 
species a, AI{Sa) is the normalization constant, g{Sbz) 
is the standard Clebsch-Gordan coefficient [l^ 
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( Q \ — / Q Q Q ■ Q Q \Q Q Q Q \ — ^ + Sb — Sy.{Sa — Sb + S)\{Sb Sg + 5^)! 

(ba + bb + O + L)\ 

X [{Sa + s,- SbzY-iSa -s, + Sbz)KSb + SbzY-iSb - SbzV-is + s.y.is - s,y.]^ 
X Eki-^)''[kKSa + Sb-s- ky.{Sa -s, + Sbz - ky.iSb + Sbz - ky.{s -Sb + s,- Sbz + ky.{s -Sa~ Sbz + ky]-\ 

(13) 

where k is an integer such that the arguments in the factorials are non-negative. The summation over k reduces 
to only one term in the maximally polarized ground states in each parameter regime, and in the ground states for 
c"-^ > 0. 

li Sz — S, we have 



diSbz) — {Sa^S Sbzi Sb, Sbz\Sa, Sb, S, S) 



= (-1) 



(25 + msa + Sb~ sy{Sa - Sbz + sy{Sb + Sbzy 



iSa + Sb + s + iy.{Sa -Sb + sy.{Sa -s + SbzyiSb - Sbzy 



(14) 



which gives giv{Sbz) in R-ef. [IjI by substituting Sa, Sb and S as Na, Nb and n, respectively. 
If 5 = 5° - 5'' > 0, we have 

9 {Sbz) — {Sa, Sz Sbz', Sb, Sbz\Sa, Sb, Sa Sb, Sz) 



(-1) 



Sbz 



{2Sa - 2Sb + iy.{2Sby{Sa + Sz- Sbzy{Sa -Sz+ Sbzy 



(15) 



i2Sa + msb + Sbzy{Sb - Sbzy{Sa -Sb + Szy{Sa -Sb- Szy 

In particular, for Sz = S = Sa — Sb, 

\{2Sa -Sb + l){2Sby{2Sa -Sb- Sbzy 



g{Sbz) — {Sa, Sa — Sb — Sbz', Sb, Sbz\Sa, Sb, Sa — Sb, Sa — Sb) = (— 1) 



li S = Sb - Sa>0, we have 

g{Sbz) — {Sa, Sz — Sbz', Sb, Sbz\Sa, Sb, Sb — Sg, Sz) 



{2Sa + msb-Sbzy 



(-1) 



Sbz 



{2Sb - 2Sa + iy.{2Say.{Sb + Sbzy-{Sb - Sbzy 



{2Sb + mSa +Sz- Sbzy.iSa -Sz+ Sbz)KSb - Sa + 5,)!(5b -Sa- 5.)! 

In particular, for Sz = S = Sb — Sa, 

' {2Sb - 2Sa + l)i2Say.{Sb + Sbzy 



9{Sbz) — {Sa, Sa — Sb — Sbz', Sb, Sbz\Sa, Sb, Sa — Sb, Sa — Sb) — (" l)"^*"' 

I 



{2Sb + l){2Sa - Sb + Sbzy 



(16) 



(17) 



(18) 



Furthermore, we can know the structure of 
\Sa,Sb,S,S) in units of singlet pairs of atoms, us- 



J 



ing the discussions in Ref. jlf 
U Na + Nb - S is even, 



ic c c c\ _ ^ ^ fQi.o.i , tQo,i,i „ f Qi,i,o„-|- Q2,o.o„-f Qo,2,0| , 

\ba,bb,b,b) — 2_^A[l^l^lfi, 1^2,0,0, ^0,2,0)O'i 0-^ ^1,1,0 ^2,0,0 ^0,2,0 1^/ 



where 



0li,o = al^-i-«X + «-i^ 



(19) 



is the creation operator for a interspecies singlet pair 62 o o = ^^d 62 = ®& ^'^^ creation operators for 
intraspecies singlet pairs, Qi,o,i and Qo,i,i are given by 



Qi,o,- 



Q2,0fi — Qo,2,0, 
+ Qo,2,0 ~ Q2,0,0, 



(20) 
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A special case is the singlet state under Na = Nf,, for which it is clear that Qi^.i = Qo,i.i = 0, that is in each term 
in the superposition, all atoms are paired up. 
UNa + Nb~S is odd, 



ic c c c\ _ /o ^ ^ t ^Qi.o,i^Qo.i,i + i if Qi,i,o | Q2,o.q -f Qo.2,0, , 

\i>a, Ob, J, j) — ^Wi,1,0j <^2,0,0, VO,2,OJao'^l "1 "14,0 *^2,0,0 *^0,2,0 1^/ 

/t'/'/o /n /n ^ t'^^-°'^^"'"i,tj,t'^°'^'^ot '^^•^'"ot Q2,o,o„-(- Qo,2,0| , 

+ ^ (,^1,1,0! V2,0,0i VO,2,oj«l Oq"! ^1,1,0 "2,0,0 "0,2,0 1^/' 



where (5i,o,i and (5o,i,i are given by 



(21) 



n.0,l — 2 ^2,0,0 — V0,2,0, /r,r,N 

^0.1,1 — 9 r V0,2.0 ~ V2,0.0- 



In both and (PT|) . the summations are over Qi^.o, Q2,o.o and Qo,2.o- The coefficients A and ^' are determined 
by constraints 

^a\^a,T Sb, S, S) ~ S'Q(S'a + 1)|S'q, 5^, 5*, 5*), (O'W 

Sl\S,,Sb,S,S) = Sb{Sb + \)\Sa,Sb,S,S). ^^^^ 



where 5^ = (ajai — aLift-i)^ + 2(aQ)^Q;iQ;_i + 

2a|Q!Liao + 2a|!)aoQ!ia| + 2Q;Jaoa_iQ;l]^. 

For a ground state which is a direct product of the 
states of individual species, one can resort to the com- 
posite structure of a single species of spin-1 atoms |5j, 
and does not need to use the result here. For an entan- 
gled state of the form of l^^^, , \S^ - S^\, \S^-S'^\), 
one can use the above result. 

Now consider a particular ground state \Na,Nb,Na — 
Nb, Na — Nb) ■ As now Na + Nb — S is even, we know that 
\Na,Nb,Na~ Nb,Na~ Nb) IS givcu by Eq. (HH), with 

I Ql,0,l = Na — Nb + (52,0,0 ^ Qo,2,Oj /2^\ 
\ Qo,l,l = Qo,2,0 ^ Q2,0,0- 

V. SUMMARY 

Under the single spatial mode approximation for each 
atom, the many-body Hamiltonian of a mixture of two 
species, labeled as a and 6, of spin-1 atoms with both in- 
traspecies and interspecies spin exchanges can be reduced 
to that of two coupled giant spins, with the intraspecies 
and interspecies spin coupling strengths c°, and c"'', 
as given in ((!]) . As listed in Table HI we have determined 
the ground states in nearly all parameter regimes, ex- 
cept a small patch where there are too many possibili- 
ties, by minimizing the energy under various conditions. 
These states can be rewritten in terms of bosonic de- 
gree of freedom, as discussed in Sec. |lVl The states can 
always be written as superpositions of various configura- 
tions of intraspecies and interspecies spin singlet pairs, 
together with some individual atoms with z-component 
spin 1, the total number of which is 5*2. For a many-body 
singlet ground state, all the atoms are paired up in each 
configuration. 

The many-body singlet ground states in two regimes 
deserve particular attention. One is the non-degenerate 



singlet ground state \Nb, Nb, 0,0) in the regime c*^^ > 0, 

< and <c'' < c"^ (excluding c°- = = c°-^) 

plus < and c"^ < c" < ^^^^Trrf-- This singlet state 
exists even in the generic case oi'Na > Nb- In the special 
case of Na = Nb, its regime only expands to < c°'^ and 

< d'^ plus < 0, and d^ < < ^^^j^t+r' P^^^ 

cb < m^±nf. while < 0. 

The other regime of many-body singlet states is c° = 
= c"'' > 0, where the ground state can be any 
15^,5'^, 0,0) with a legitimate value of S^. Therefore 
in the generic case of Na > Nb, in the regime c"'' > 0, 

< c"'' and < < c"'' (including c" = c'' = c"^) 

plus c'' < and < c" < ^^^^^l^j— , we always have 
S = Sz = 0, even though the magnetic field is absolutely 
zero. This robustness is useful for experiments. 

At c"' — c'' = c"'', there is a continuous quan- 
tum phase transition from the global singlet state 
\Nb, Nb, 0,0) to the product of two single-species sin- 
glet states |0, 0)a|0, 0)b, which is the ground state in the 
regime d > 0, c'' > and c^c'' > (c"'')^. 

Focusing on the maximally polarized states with 5*2 — 
5™, as can be picked out by an infinitesimal magnetic 
field, we have described the phase-diagrams on c° — c'' 
planes for zero, negative and positive c"'', respectively, 
thus obtaining the three-dimensional phase diagram in 
c° — c'' — parameter space. 

For = (FIG. [3]), the four quadrants correspond to 
the products of two antiferromagnetic states |0, 0)a|0, 0)b, 
of a ferromagnetic state and an antiferromagnetic state 
\Na,Na)a\0,0)b and \0,0)a\Nb,Nb)b, and of two ferro- 
magnetic states \Na, Na)a\Nb, Nb)b, respectively. Each 
ground state is a distinct quantum phase, with the 
and axes as the phase boundaries, on which there are 
degenerate ground states, two of which belong to the two 
bordering phases, respectively. 

Another kind of phase boundaries appear in the case of 
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< (FIG.[T]). With the increase of |c°^|, there grows a 
surface c'' = in the part > and a surface c° = 

in the part c > 0. These two surfaces, together 



2N, 



b£_ 



and 



with c — c plane, (c" = 0) in the part c > -2n~TT' 
c°- - c°-^ plane (c'' = 0) in the part > |^^, act as the 
phase boundary separating |0,0)a|0,0)b from the other 
phase, to which the phase boundary itself also belongs. 

The latter phase for c"** < consists of five ground 
states, which are continuously connected on four bound- 



aries, c 



-2Afi,c'" 

2Na + l 



and 



and 



-2NbC°-'' for Cb < 0, and 

-2N;:pr '^''^ ^ - -2A^ac"'' for < 0. These 
four boundaries of continuous connection converge re- 
spectively to negative halves of c" and axes as c""^ — >■ 0. 
Especially, for a given c"'' < 0, the crossover regime 
of \Na, Na)a\ni,Tii)b disappears as c"'' — ^ 0. Conse- 
quently, \Na,Na)a\Nb,Nb)b and |iV^,iV„)a|0,0)6, which 
can be continuously related through \Na, Na)a\ni,ni)b, 
become discontinuous neighbors as c""^ — ^ 0. Simi- 
larly, the crossover regime of \n2,n2)a\Nb, Nh)b also dis- 
appears as c"'' — > 0. Consequently, \Na, Na)a\Nb, Nb)b 
and \0,0)a\Nb, Nb)b, which can be continuously related 
through \n2,n2)a\Nb, Nb)b, become discontinuous neigh- 
bors as c"''' — 0. This conforms with the c° — c'' phase 
diagram for c"'' — 0, in which the negative c° axis and 
the negative axis are phase boundaries. 

With the increase of c""^ on the positive direction, 
there also grow a number of boundaries where the 
ground states are continuously connected, consequently 
the ground states in the second, third and fourth oc- 
tants, plus the part of first octant with c" < and 



< c (excluding c° 



c ) are all continuously 



connected to be a single quantum phase (FIG For 
c" < 0, there are two boundaries of continuous con- 



nection c' 



b _ 2(Af„ + l)c° 
2Nb + l 



and 



2{Na + l)c'^^ be- 



tween which is the crossover regime of ground states 
\Na,n4, Na — n4,Na — rii), which is continuously con- 
nected downwards with \Na, Nb, Na — Nb, Na — Nb) and 
upwards to | A^q, A^a)a|0, 0)b. In case Na > Nb, there are 

two boundaries of continuous connection — 2n'''+i 

2N. 



and c'^ — 2n''+i ^ — ' '^^ich disappear in case 
Na = Nb. fhese two boundaries define the crossover 
regime of | n2 , A^f, , ^2 — A^f, , n2 — A*";,) , which is continuously 
connected leftwards with \Na, Nb, Na — Nb, Na — Nb) , and 
rightwards with | A^b, A^b, 0, 0). For c'' < 0, there are two 



boundaries of continuous connection c"^ 



2(Aft, + l)c° 
2Nb + l 



and 



2(A'b + l)c°'^, between which is the crossover regime of 
|n3, A'b, A'b — n3,Nb — n^), which is continuously con- 
nected leftwards with jA'b, Ab, 0, 0), and rightwards with 
|0, 0)a|A^b, A^b)fc- When c"'' — )■ 0, the crossover regimes 
all tend to disappear, therefore \Na, Na)a\0,0)b and 
|0, 0)a|A^b, A'^b)?) tend to occupy the second and fourth 



quadrant of the 



cross-section, as for c' 



,ab 



0. 



But unlike the case of c"'' < 0, there is a discontinuity in 
the fourth quadrant between | A^q, A^b, A^a — Nb, Na — Nb) 
for c"'' > and \Na, Na)a\Nb, Nb)b for = 0. In the 



first octant, with the increase of c°^, the boundary of 
|0, 0)Q|0,0)b grows from the positive c° and axes to 
c"c^ = (c°^)2, c° > and c'' > 0. 

Interspecies spin coupling c"'' provides an extra pa- 
rameter controlling quantum phase transitions. We have 
given two examples. First, the discontinuity between 
ground states for c"'' > and c"'' < when c° < and 

< implies a first order quantum phase transition in 
varying the sign of c"''. Second, by tuning c"'' from zero to 
negative and then back to zero, the first order quantum 
phase transition between negative and positive values of 
an intraspecies spin coupling can be circumvented. 

With the experimental experiences on single-species 
spinor gases [g, and on scalar mixtures [§], where in- 
terspecies spin exchange has been observed as a disad- 
vantage, and with the development of heteronuclear Fes- 
hbach resonance [l9j . future experimental realization of 
a spinor mixture with interspecies spin exchange, with 
the ground states described here, is expected. 
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Appendix A: 5^,5^ and S™ in the case c"' < 



In this appendix, we find S"™, S"™ and S"\ in which E 
is minimal, in parameter regimes with c"'' < and in the 
absence of a magnetic field. In the discussion, E always 
represents the energy as low as can be determined in the 
regime under discussion, i.e. the meaning of E keeps 
updating. 

With c"'' < 0, i? is minimal when S = Sa + Sb- Hence 

E = ^SaiSa + 1) + jSbiSb + 1) + c'^'SaSb. (Al) 



Thus 



dE 
dS'a 



c'^Sa + c'^^Sb 



2 ' 



dSb^ 2- 



(A2) 



(A3) 



Several different cases are considered in the following. 



1. < 0, c*" < 



In this regime, of which regime I in Ref. [16| is a subset, 
^ < 0, ^ < 0, hence S'^ = A^„ = A^b,^'" - 

A^a°+ Nb. 
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FIG. 6: Possible values of Sa and St are within the rectangu- 
lar defined by < Sa < Na and < St < Nb on Sa-St plane. 
The dashed line is the stationary line 4§- = in the case 



c"" < and < c" < 



2JV„c°° 
' 2JVi, + l 



The solid line is the station- 



ary line in the case c" < and 



< c" < -2iV„c'" 



c° < 0, c*" > 



In this regime, < 0, hence 



where Int(a;) represents the integer closest to x and in the 
legitimate range of St] here < ni < Ni,. <S"" = Na + ni. 

c. > -2NaC'''' 



If > —2NaC'^^, all points {Sa,Sb) in the rectangular 
satisfy ^ > 0. Therefore S*™ = 0, S*™ = S"" = Na- 
The same solution is obtained if = —2NaC°'^. 



3. c° > 0, c*" < 

One simply exchanges the subscripts or superscripts a 
and b in the preceding subcase. Thus there are also three 
possibilities. 



a. < c° < 



Na+1/2 



= Na, SJ^ = Nb, 5" ^ Na + Nb. This regime can 
also be combined with lA 11 with the same solution. 



E{Na, Sb) = -SbiSb + 1) + c"'''NaSb + const. (A4) 

We represent all the values of Sa and Sb as points 
{Sa,Sb) within the rectangular defined hy < Sa < Na 
and < S'b < iVfc on Sa-Sb plane (FIG. E)- ^ = de- 
fines a stationary line. The points above this line satisfy 
> 0, while the points below it satisfy < 0. One 
can see three possibilities. 



b. < < -2iVi,c"\ 



iVa + l/2 



5*™ = n2, with 



Ic"''! 1 
n,^lnt[^-ANb^-l 



(A6) 



< c" < 



2JVoe°'' 
' 2JV6 + 1 



with < 712 < Na, Sir = Nb, 5™ = na Nb 



The stationary line, depicted as the dashed line, crosses 
with the line Sb = Nb- Hence all points with Sa — Na 
satisfy If < 0. Consequently 5™ = Na, 5^ = Nb, 
gm _ Note that this regime can be combined 

with case I A 11 where the solution is the same. 



The stationary line, depicted as the solid line in FIG. El 
crosses with the line Sa = Na- The crossing point gives 
the minimal energy. Hence 5™ = Na, S'™ = ni, with 

lc"''l 1 
ni=Int[^A^,--], (A5) 



4. c" > 0, < c*" < -2iVac"\ c"^c'' > {c"^^f 

First we consider < c'' < — i^^-r-r. Then the sta- 

tionary line — crosses with the line Sb = Nb, as 
shown as the dashed line in FIG [SI The crossing points 
with Sb = and Sb = Nb divide the rectangular to three 
regions, for which the minima of E are 
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E 



^Sa{Sa + 1), 
2c'' '^a 

^SaiSa + 1) + C-''NbSa + ^NbiNb + 1), 



if < < 
if - 
if - 



< 5„ < - 



(jVt + l/2)c'' 



(A7) 



r 



5*6 in the first interval of Sa- In tfie second interval, 
and Sb satisfy dE/dSb = 0. In the third interval, 
Sb = N},. In each of these three intervals, E increases as 
Sa increases. Therefore, = S*^" = 5™ = 0. 



Now we consider 



2Na 



2Nt + l 



< -2iVac"^ Then the 



stationary line = crosses with the line Sa — Na, 



as shown as the solid line in FIG. El The crossing points 
with Sb — and Sa — Na divide the rectangular to two 
regions. The minima of E are given by 



E : 



"-SaiSa + 1), 

D„ - 



2c'' 



-Sa 



if 0<Sa< 



2cai> ' 



if 



20" 



<Sa< Na. 



(A8) 



S'b = in the first interval of Sa- In the second region, 
Sa and Sb satisfy dE/dSb — 0. Again it is found in each 
region, E increases as Sa increases, therefore 5"™ — S™ — 
S"" = 0. 

Therefore, throughout this parameter regime, we have 
S^ = Sj;' = S"" = 0. 



c' > 0, < < -27V6c''\ c'^c' > {c''''f 



Exchanging labels a and b in the preceding subcase, 
we obtain 5*™ = 5^™ = 0. 



6. > 0, c' > -2Nac'''' 

All points satisfy dE/dSb > 0. Hence SI'' = 0. Conse- 
quently, one also obtains, from Eq. (|A1[) . S*™ — 0. 



9. c'' > 0, < c" < -27V6c"\ cV* < (c"'')2 

Exchanging labels a and & in the preceding subcase, 
we also obtain 5™ = 5^" = fi-™ = 0. 



10. 



< < -1^, < c'' < 



2JVaC°'' 

" 2iVi, + l ' 



C C < (C J 



-B is as given in (IA7|) . In the first interval of Sa, the 
minimum of E is at Sa — 0. In the second interval, 
the minimum is either at Sa = — ^'^''^J^C^^'^ , which also 
belongs to the first interval, or at Sa — ~2c^' '^hi'^f^ 
also belongs to the third region. From c" < — ^^^1/2 ' 
have 71,2 > Na- Consequently, the minimum in the third 
interval is at Sa = Na, moreover, it is straightforward to 
see E{Sa = Na) < E{Q). Therefore, S^' = Na, = Nb, 

5™ ^Na+ Nb. 



7. c' > 0, c" > ~2Ntc°- 



Exchanging labels a and b in the preceding subcase, 
we have S"™ = 5^ = 



Appendix B: 



> 



In this appendix, we find out 5™, 5™ and 5"", in which 
E is minimal, in parameter regimes with c"'' > and m 
absence of the magnetic field. Without loss of generality, 
assume Na> Nb. 



Now E is as given in ([M]) . With c°c'' < (c°^)^ 
dE/dSa > in the the second region. Consequently 



S" 



sj: 



0. 



In this case, Sa and Sb both disappear in E. Hence E is 
minimal when 5 = 0. Thus 5*™ = S*™ for any legitimate 



value of 5^™, i.e. S*™ = 5^ = 0, 



,Nb 
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Sb' 










Nb 




(A^A, Nb) 






I /'' 


/^''''' 




■ 






II 






O/ 




N 


a 


Sa 



Therefore, the minimum of E in region II hes on the 
border line Sa — St for < Sa < Ni,, and on the line 
Sb = Nb for Nb<Sa< Na. 

Therefore, the minimum of E in the whole rectangular 
must lie on the line interval Sb — Nb, Nb < Sa < Na, on 
which 



E = y5a(5a + 1) - c'^'NbSa + COUSt. 

There are three possibilities. 



(B6) 



FIG. 7: Possible values of Sa and St are within the rectangu- 
lar defined by < Sa < A''a and < St < Nb on Sa-Sb plane. 
The solid line represents Sa ~ Sb, above which is region I 
while below which is region II. The dashed line is the station- 
ary line = in region I in the case c" > c"'' > c*" > 0. 
The broken line is the stationary line in region II in the case 
> c"'' > c" > 0. 



In any other cases, as considered in the following, E is 
minimal when 5 = 1 — S'fc I . 



This includes the range that c" < 0. S^" = Na, SJ^' 
Nb, S"' ^Na- Nb. 



2Na+l 



2JVi, + l 



We have 5™ = na, S"^" = ^6, 5™ = ^2 - iVb. 



iN^c-^" ^ „a ^ ab 
2iVb + l — ^ — 



2. c" < c"", c" < c"", except c" = c" = c" 



We have SJ! 



sr 



Nb, S" 



As shown in FIG. [3 we divide the whole region of Sa 
and ^f, to two regions by the line Sb — Sa- In region I, 
E is minimal when S = Sb — Sa, i.e. 



E = ^Sa{Sa + l) + ^SbiSb + l)-C^'SaSb-c"'''Sa. (Bl) 

Hence 

— = c'^Sa - c'^'Sb y - C"^" < 0. (B2) 

Therefore, the minimum of E in region I lies on the bor- 
der line Sb = Sa, on which 



\ Sa{Sa + l) 



(B3) 



3. > c"'' > c'' > 0, c^c'' > {c''''f 

In region II, < 0, as shown in Eq. (|B5p . Hence 
the minimum of -B lies on the border line Sa = Sb for 
< S'a < Nb, and on the line Sb = Nb for Nb < Sa < Na- 
When Sb = Nb, E is given by (|B6I) . for which we now 
have 

— = c-^Sa - c-'Nb + y > (c-^ - c''^)iVfc + y > 0- (B7) 

Therefore in the line interval Sb = Nb for Nb < Sa < Na, 
the minimum of E locates at {Nb, Nb). Consequently the 
minimum of E in the whole rectangular lies in region I, 
in which E is given in (|Bip . thus 



which reaches its minimum at {Nb, Nb), as c°' + c''—2c'^'' < 
0. 

In region II, S = Sa — Sb, 



E = -Sa{Sa + l) + jSb{Sb+l)-C-''SaSb~c''''Sb. (B4) 



Hence 



^ = c'Sb - c'^'Sa -f ^ - c-^^ < 0. (B5) 

OJb ^ 



dSb~ ' + 2 



(B8) 



= defines a stationary line which crosses with the 
lines Sa = ^6 and Sb = Nb at (^(jsi-rpry , 2(c°fc-c'') ) '^^^ 
^£_(7Vb_fi/2) ^ ^^-j j-ggpgctively, depicted as the dashed line 
in FIG. El 

If ^ < ' t'^^ crossing point with Sb = Nb lies 

in region I, then the minima of E can be found to be 
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E 



c''+c°-2e' 



2 Sa{Sa + 1), 

-ab-.2 „ a ab 

q2 I c — c c 



if < < 
if 



2(cai.-c'') 



< < 



ab )\L e „ab Q , f AT. i 1 \ ;f c''{Nt + l/2) ^ g < Nh 



c''(JVt, + l/2) 



^SaiSa + 1) - C-'N^Sa ~ C^'Sa + ^Nh{Nt + 1), if 



(B9) 



Thus 



dE 
dS~a 



(c°+c''-20(5, + i), if 0<Sa< 



2(c"''-c'') 



-5a 



if 



<5a < 



c''(JVt, + l/2) 



„ab 



if '^^^^<S^<N, 



(BIO) 



As c^c^ > (c°'')^ c" + c*" > 2VJ^ > 2c°^ we find that 
> in all the three intervals. Consequently — 
sf = S*™ = 0. 

If ^ > -jq^jp^2- §§l > in the whole region 

II. Then the minima in the whole region I are given by 
the first expression in (|B9[) . One again obtains = 

grn ^ gm ^ q 

4. C*" > c"' > c" > 0, c^c'' > (c"'')2 

One still uses FIG. [7l In region I, > 0, hence the 
minimum of E lies on the border line 5*^ — Sb- Thus the 



minimum of E in the whole rectangular lies in region II, 
in which E is given in (jB4|) . and thus 

— ^c-Sa-c'''Sb + -. (Bll) 

= defines a stationary line, shown as the broken 
line in FIG. [71 which crosses with the line Sa — Sb at 
( £! C ^ 

^ 2(c"'>-c'') ' 2(c"'>-c") ^ 

If ^ ^ jy +1/2 ' thfi — also crosses with the 

boundary Sb — Nb . Then the minima of E can be found 
to be 



^"^2 — Sa{Sa + 1), 



iSaiSa + 1) - C^^'NbSa + ^ Nb{Nb + 1) - c'^'^Nb, 



if < S'a < 
2f^-l), if 



2(cai>_c») 



if s^iiVb - i < s-, < iv. 



(B12) 



> (c"'')2, it is found that ^ > in all the 



With c^-d 

three intervals. Consequently 5" 



SI! 



S'' 



0. 



If 



< 



Nb 



the stationary line = crosses 
Na- Then the minima of E in 



Na + l/2 ' 

with the boundary Sa 
region II is given by the first two expressions in (jB12l) . 
with the second interval of Sa extended to Na- One ob- 



tains the same result ST 



0. 



6. c" > c"^ >0, c'' <0 



The reasoning given in case IB 31 that the minimum of 
E must lie in region I is still valid. Using (jBSp . now 
dE/dSb < 0, as c'' < 0. Therefore S^ = Nb, for which 



E = y S„(5a + 1) - c'^'NbSa - c'^'Sa + const, (B13) 



with < Sa < Nb- There are three possibilities. 



This regime has been discussed in [l&l, with result 



sj; 



S" 



0. 



The cases IB 31 IB 41 and the present one can be written 
altogether as c°c'' > {c"'^)'^, with the same result 5™ = 



SJ! 



0. 



a. e°''<c''< ^^%+'2r , c" < 

— JA/f, + l — 



S™ = N'', SI"" = Nb, 5™ = 0. 
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b. ^^^^+^\"'"' < < 2{Ni, + l)c'^^ c*- < 
5™ = ns, where 

"3 = Int[ -], 

with < 713 < Nb. SY" = Nb, S"" ^Nb- ns- 



(B14) 



2(Na + l)c' 



2Nh 



^ < c" < 2(Ar<, + 1)0"% c" < 



5,"^ = iV„, 5," = 714, where 

"4 = Int[ -f, -] 

with < n4 < Nb. S"' = Na- n^. 



(B16) 



= 0, Sr = S'' 



c. > 2{Nb + l)c"\ c'' <0 
Nb. 



c. c'' > 2{Na + l)c'"\ c" < 
S"^ = 0, S"^ = 5"^ = Na. 



7. J' > c"*" > 0, < 



As c° < c"'' , the reasoning in case IB 41 that the min- 
imum of E must he in region II is stih vahd. Then 
dE/dSa < 0. Therefore ^ Na, for which 

E = —Sb{Sb + 1) - c'^^NaSb - c^^Sb + const. (B15) 
There are three possibihties. 

a. e'"'<c''< ^'^°^+;f\ c°<0 
5^ = iV^ S'^ = Nb, 5™ = iVa - Nb. 



8. > 0, c' > 0, c" + c'' < 2c"'' 

The case of c° > 0, c'' > and c^c'' < (c"'')^ has too 
many possibihties to be calculated in detail here. We only 
consider the sub-regime c° > 0, > 0, c° -I- c** < 2c''''. 
According to Equations (jB9p and (jB12l) . there is always 

an interval of Sa with E = ^^-^A^^Sa{Sa + 1) and 
Sa = Sb. As c°- + < 20°-^, it is certain that in this 
regime, E is minimal at the upper end of this interval, 
although it remained undetermined whether at this point 
E is the minimal in the whole region of Sa and Sb. 

Hence it is determined that S*™ ^ 0, 5™ 7^ 0, conse- 
quently the ground state | S*™ , S*™ , 1 5™ - S*^"^ | , | S*^ - 5^™ | ) 
is always entangled. 
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